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Abstract 



(<~^ ' The paper is concerned with open quantum systems whose Heisenberg dynamics are described by quantum 

^SJ , stochastic differential equations driven by external boson fields. The system-field coupling operators are assumed 
to be quadratic polynomials of the system observables, with the latter satisfying canonical commutation relations. 

, In combination with a cubic system Hamiltonian, this leads to a class of quasilinear quantum stochastic systems 

^) ' which retain algebraic closedness in the evolution of mixed moments of the observables. Although such a 

, system is nonlinear and its quantum state is no longer Gaussian, the dynamics of the moments of any order are 

' amenable to exact analysis, including the computation of their steady-state values. In particular, a generalized 

, criterion is developed for quadratic stability of the quasilinear systems. The results of the paper are applicable 

,— — i' to the generation of non-Gaussian quantum states with manageable moments and an optimal design of linear 

, quantum controllers for quasilinear quantum plants. 

C : I. INTRODUCTION 

a . 

Q-i' The paper is concerned with open quantum systems whose Heisenberg dynamics are described by 
quantum stochastic differential equations (QSDEs) driven by external boson fields [13]. This class of 
y—i dynamical systems, whose variables are noncommutative operators on a Hilbert space, is a common 
^ source of models in quantum optics [2] and quantum control [1], [4], [5]. Central to this approach is 
\Q a Markovian description of interaction between the quantum-mechanical system and its environment, 
, whose "energetics" is specified by the system Hamiltonian and system-field coupling operators. The 
^ case where the Hamiltonian is quadratic and the coupling operators are linear in the system observables, 
^ satisfying canonical commutation relations (CCRs) [11], corresponds to an open quantum harmonic 
psj oscillator. Its importance for linear quantum stochastic control [1], [12], [15], [16] is explained by 
^ ! tractability of the dynamics of moments of the system observables, which is closely related to the 
>• ■ invariance of the class of Gaussian quantum states [14] of the system subject to external fields in the 
vacuum state [13]. This allows the linear quantum control to inherit at least some features of the classical 
^ control schemes [7], including practical computability of quadratic costs. In the present paper, we 
consider more complicated system-field interactions, described by coupling operators which are quadratic 
polynomials of the system observables. In combination with a cubic system Hamiltonian, this leads to a 
novel class of quasilinear quantum stochastic systems which extend the linear quantum dynamics above 
and yet retain algebraic closedness in the evolution of mixed moments of the observables, similarly to 
their classical counterparts in [17]. Although such a system is nonlinear and its quantum state is no 
longer Gaussian, the dynamics of the moments of any order are amenable to exact analysis, which 
includes the computation of their steady-state values. In particular, this allows a generalized criterion to 
be developed for quadratic stability of the quasilinear systems. The results of the paper are applicable 
to the generation of non-Gaussian quantum states [18], [19] with manageable moment dynamics. They 
can also be used for an optimal design of linear quantum controllers for quasilinear quantum plants. 
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II. OPEN QUANTUM SYSTEMS 

We will briefly review the underlying class of models for open quantum systems which is based on 
QSDEs. Suppose W{t) :— {Wk{t))i^ki^,„ is an m-dimensional quantum Wiener process of time ? ^ on 
a boson Fock space ^ [13], with the quantum Ito table 

dWdW^ := idWjdWk)i^j,k^„, = adt, (1) 

where the time argument of W is omitted for the sake of brevity, and Q. := (ci);fc)i<;,A:s;m is a constant 
complex positive semi-definite Hermitian matrix of order m. Here, a complex number (p is identified 
with the operator cpl,^, where denotes the identity operator on Also, vectors are organised as 
columns unless indicated otherwise, and the transpose (■)^ acts on matrices with operator- valued entries 
as if the latter were scalars. The entrywise real part of the matrix Q, denoted by 

V := ivjk)i^j,k^m ■= Ref^, (2) 

is a positive semi-definite symmetric matrix. The entries Wi{t), . . . ,Wm{t) of the vector W{t) are self- 
adjoint operators on associated with the annihilation and creation operator processes of extemal 
boson fields. The imaginary part of the quantum Ito matrix t2 in ([T]) describes the CCRs between the 
fields in the sense that 

[dW,dW^]:={[dWj,dWk])i^j^k^^ = iJdt, 7:=2M, (3) 

where [A, 5] := AB — BA is the commutator of operators which applies entrywise, and i := is 
the imaginary unit. We consider an open quantum system (further referred to as the plant) whose 
state variables Xi(0), . . . ,X„(0) at time ? = are self-adjoint operators on a complex separable Hilbert 
space J^. Any such operator ^ can be identified with its ampliation on the tensor product 

space J^0.^. Suppose the plant interacts only with the extemal fields, so that the density operator 
p{t), which describes the quantum state [11] of the plant-field (as a closed system, isolated from the 
environment) on ^ at time t, evolves in the Schrodinger picture of quantum dynamics as 

p(0 = f/(OP(o)f/(0^- (4) 

Here, U (t) is a unitary operator on the Hilbert space J^f ® initialized at the identity operator 
U{0) =Ijf(g,,^, and (■)^ is the operator adjoint. The initial quantum state of the plant-field composite 
system is assumed to be the tensor product 

p(0) := 07(0)01; (5) 

of the initial plant state 07(0) on and the pure state v := |0)(0| of the extemal field associated 
with the vacuum vector |0) in where use has been made of the Dirac bra-ket notation [11]. In the 
Heisenberg picture, an observable i^(?) on ^ ®^ evolves in a dual unitary fashion to dH): 

^(0 = t/(0^^(o)f/W, (6) 

with the duality being understood in the sense of the equivalence 

E^(0:=Tr(p(0)^(0)=Tr(p(0^(0)) (7) 

between two representations of the quantum expectation. The unitary operator JJ (t) itself is driven by 
the intemal dynamics of the plant (which the plant would have in isolation from the surroundings) and 
by the plant-field interaction. In the weak interaction limit, which neglects the influence of the plant on 
the Markov stmcture of the field in the vacuum state, a wide class of open quantum systems is captured 
by the following QSDE: 

dU = -{{iH + h^ah/2)dt + ih^dW)U. (8) 
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Here, H is the plant Hamiltonian, and h := {hj)i^j^m is a vector of plant-field coupling operators, with 
hj pertaining to the interaction between the plant and the jth external field. Both H and hi,...,hm 
are self-adjoint operators on which are usually functions of the plant observables Xi(0), . . . ,X„(0). 
The term h^AW = ^Li^k^Wk in which can be interpreted as an incremental perturbation to the 
plant Hamiltonian H due to the interaction with the external fields, is an alternative form of the more 
traditional representation 

h^dW = i{L^d£/*-L^ 

through the m/2-dimensional field annihilation and creation operator processes £/{t) 
Fock space ^ with the quantum Ito table 



and 



(9) 

)* on the 



4,/2 





dL 



(10) 



so that \ds^,ds^^\ := ds^ ds^^ - ds^^ djj"^ = /„/2dL Here, m is assumed to be even, (•)''' := 
denotes the transpose of the entrywise adjoint (■)*, and /,■ is the identity matrix of order r. In application 
to ordinary matrices, (■)^ reduces to the complex conjugate transpose (•)* := ((■))^- The representation 
Q corresponds to the case when the scattering matrix [13] is the identity matrix. The vector L 
{Lk)i<^k^m/2 consists of linear operators on J^, which are not necessarily self-adjoint. The relation of 
h, W with L, £/ is described by 



1 
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where ® denotes the Kronecker product of matrices, so that the corresponding quantum Ito matrix Q in 



© is fl = 
form h^Q.h 



Im + i 



'm/2 



/2. In general, the operator h^Q.h := Jl'fk=i ^jkhjhk, which takes the 



in the case ([I0l)-([l2l), is self-adjoint since {h^^hy =Z'li^^iaJj^hkhj 



h^a*h = h^Qh. The formulation using h and W, in principle, allows W to be an additive mixture of 
the quantum noise with a classical random component such as the standard Wiener process. In view of 
([3]), the "magnitude" of ImQ. = 7/2 in comparison with ReQ. (measured, for example, by the spectral 
radius r(ImQ(Ref2)^^) which is well defined and is strictly less than one if Q 0) indicates the relative 
amount of "quantumness" in W. This setting reduces to the classical noise situation if the matrix 
in (dl) is real, in which case U (t) becomes a random process with values among unitary operators on 
see [6] and [13, pp. 258-260]. The general situation is treated by applying the quantum Ito 
rule d(77Q = (dT])^-|-77d^ 4- (dr])d^ and using ([8]) along with the unitarity of lJ{t) and commutativity 
between the forward increment dW and the adapted processes. This yields the following QSDE for the 
density operator p{t) in dH): 

-/([//, p\dt + [/i^, p]dW) + \^{Q}C{p))dt, 



dp 



(13) 



which is referred to as the stochastic quantum master equation [2] in the Schrodinger picture. Here, 
use is made of a self-adjoint operator tr(Q^C(p)) :=LJyt=i '^jkCjkip) on ^0^, with tr(-) denoting 
a "symbolic" trace (to be distinguished from the complex- valued trace Tr(-) of an operator), where the 
matrix C(p) := (Cjyt(p))i^j,/t<ra has operator-valued entries 



Cjkip) := hkphj - {hjhkp+phjhk)/2. 

Note that Cjk{p)^ = Ckj{p) for any self-adjoint operator p in view of self-adjointness of h\. 
this, together with Q.* = Q., ensures the self-adjointness of tr(Q^C(p)). 



(14) 

,hm, and 
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III. DECOHERENCE OPERATOR 

In the Heisenberg picture, an observable ^ (t) on the composite Hilbert space ® which 
undergoes the evolution Q, satisfies the QSDE 

d^=i{[H,Qdt+[h,Q'^dW)+^{^)dt. (15) 

Here, both the plant Hamiltonian H{t) and the vector h{t) of plant-field coupling operators are also 
evolved by the flow However, they depend on the vectorX(r) := (Xyt(f))i<;fc^n of the plant observables 
in the same way as H{0) and h{0) do on X{0). Also, =Sf denotes the Gorini-Kossakowski-Sudarshan- 
Lindblad (GKSL) superoperator defined by 

m 

^(^) ■.= tr{a'D{^)) = £ coj.Djkia (16) 

where D{£^) := {Djk{£i))i^j^k^m is a matrix with operator- valued entries 

Djki^) := hj^hk - {hjhk^ + ^hjhk)/2 

= ihj[^M + [hj.^]hk)/2. (17) 

The superoperators Dji^ are dual to Cj^ from ([14] ) in the sense that Tr(Cyyt(p)<^) = Tr(pDj^((^)). The 
superoperator matrix D acts on an observable ^ as 

D(^) = (M^,/.T] + [/,,^]/^T)/2. (18) 

The superoperators Cjk and their duals Djk, defined by (fT4)) . (flTI) . play an important role in the generators 
of quantum dynamical semigroups [3], [8]. One of such semigroups govems the evolution of the reduced 
plant density operator (n{t) obtained by "tracing out" the quantum noise in (fT3T ) over the field vacuum 
state V which yields an ODE (D = —i[H, U5] -|-tr(Q^C(Q7)). The generator of the corresponding semigroup 
in the dual Heisenberg picture is i[H, ■] + where the superoperator =Sf, given by (fT6l) . is responsible 
for decoherence [2] understood as the deviation from a unitary evolution which the plant observables 
would have alone in the absence of interaction with the environment. It is convenient to apply the QSDE 
(fT5l) entry wise to the vector X{t) of plant observables as 



dX = Fdt + GdW. (19) 

The n-dimensional drift vector F{t) and the dispersion {n x m)-matrix G{t) of this QSDE, defined by 

F :=i[H,X]+^{X), G:=-i[X,h^], (20) 

are completely specified by the plant Hamiltonian H, the quantum Ito matrix of the field process W 
from ([T]), and the vector h of plant-field coupling operators. 

Lemma 1: The GKSL superoperator (fT6l) . applied to the vector X of plant observables, can be 
computed in terms of the dispersion matrix G from (|20|) as 



1 m 

^{X) = -(GJh + i £ CDjk[hj,gk]). (21) 

Here, the matrix J is defined by ([3]), and gi, . . .,g,„ denote the columns of G: 

gk = i[hk,X]. (22) 
Proof: In view of the antisymmetry of the commutator, it follows from (l20l) that 

m 

£ 03jk[hj,X]hk = -[X,h^]ah = -iGQ-h. (23) 
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Therefore, since the quantum Ito matrix Q in ([T]) is Hermitian, then 

m m m 

( £a)j/t/v[^,%]) = Y.Wk[hk:^]^j= £ Oijk[hj,X]hk. (24) 

j,k=l j,k=l j,k=l 

By combining with it follows from (dH), ^7} that 

^{X)=i{{Gah)*-Gah)/2. (25) 
In terms of the columns of the dispersion matrix G in (|22|) . 

m 

=Gm+ Oijk[hj,gk]- (26) 

By substituting (|26l) into (|25l) and using the relationship — Q) =7 from (|3]), it follows that .^{X) — 
i(G{Q.-Q.)h + Y:jj^^^ (Ojk[hj,gk])/2 = GJh/l + iZJj^^^ (Ojk[hj,gk\/2, which establishes (ED- ■ 
In the next section, we will employ Lemma [T] in order to review the computation of the drift vector 
F and the dispersion matrix G for a class of linear open quantum systems. 

IV. LINEAR PLANT-FIELD COUPLING 

Omitting the time dependence, suppose the plant observables Xi,...,Xn, which are assembled into 
the vector X, satisfy CCRs 

[X,X'']:={[Xj,Xk]),^j,^^ = i@, (27) 

where := (0ji)i<;j,/t<n is a constant real antisymmetric matrix of order n (we denote the space of such 
matrices by A„). In this case, if the plant-field coupling operators hi,. . .,hm are polynomials of degree 
r in the plant observables, then the entries of the dispersion matrix G in (|20l ) and the vector ^(X) in 
(|2T)) are polynomials of degrees r — 1 and 2r — 1, respectively. This property follows from the reduction 
of a polynomial degree under taking the commutator with the observables due to the CCRs dTT] ): 



r 



[^k.Xi] = i ^ dkje'^ki...kj^ikj+i...kr^ (28) 

7=1 

where 

Zk:=Xk^x...xXk^ (29) 

denotes a degree r monomial of the plant observables specified by an r-index k := (fci,...,fc^) e 
{!,...,«}'', with the order of multiplication being essential in the noncommutative case. The right- 
hand side of (l28l) is a polynomial of degree r — 1 . In the plant-field interaction model which is used in 
linear quantum control [1], [5], [12], [15], the vector h of coupling operators depends linearly on X in 
the sense that 

h:=MX (30) 

for some matrix M E R*"^". In this case, the dispersion matrix G in (|20l) becomes a constant real matrix, 
since 

G = -i[X,h^] = -i[X,X^]M'^ = 0M^, (31) 

where the bilinearity of the commutator is combined with the CCRs (|27]) . In view of Lemma [H this 
allows ^{X) to inherit from h the linearity with respect to the plant observables. 
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Lemma 2: In the case of CCRs (|27l) and linear plant-field coupling (|30l ). the vector X of plant 
observables satisfies a QSDE 

dX = +KX)dt + Bm, (32) 

where the matrices K G M"^" and 5 G M"^™ are related to the the matrix J from (|3]) by 

:= 57M/2, B := 0M^. (33) 
Proof: The constancy of the dispersion matrix G, computed in (|3T1) . implies that the commutators 
on the right-hand side of (1211 vanish, that is, [hj.gk] = for all 1 ^ j,k^m. Therefore, 

^{X) = GJh /2 = ©hfjMX /2. (34) 

The QSDE ([321) can now obtained from by substituting dlB and into ([201) and using the 
notation ([33]). ■ 
Since K and 5 are constant matrices, the QSDE (|32|) may acquire nonlinearity with respect to the plant 
observables only through a nonquadratic part of the plant Hamiltonian H. Indeed, if // is a quadratic 
polynomial, that is, 

^ := E (U + ^ E rJkX^X^ = {y+RX/2fX, (35) 

where 7:= (7^)i^yt^;j ^ IR" is a given real vector, and R := {rjj^)i^j i^<^„ is a given real symmetric matrix 
of order n (we denote the space of such matrices by §„), then the commutator identities [11, Eq. (3.50) 
on p. 38]) and the CCRs ^ imply that 



m^X] = - / E (rk[X,Xk] + ^ £ rjk[X,XjXk]) 



.0y_l £ rjk{[X,Xj\Xk + Xj[X,Xk]) 

--®y+ \ E Ofc(0,X, + 0,Z,) = 0(7+i?Z), (36) 



where 0^ := {Okdi^^ki^n denotes the Ah column of the CCR matrix 0. Since the right-hand side of (|3& 
depends affinely on X, then (l32l) becomes linear with respect to the plant observables: 



dX = (AX -f 07) At + Em, A=®R + K, (37) 

which corresponds to an open quantum harmonic oscillator [1], that is, a common model employed in 
linear quantum control. 

V. ALGEBRAIC CLOSEDNESS IN MOMENT DYNAMICS 



The linearity of the QSDE (|37l) ensures algebraic closedness in the evolution of the mixed moments 
of the plant observables, defined by 

M^(0:=ES^(0 (38) 

in terms of the quantum expectation (|7]) applied to the monomials (|29|) . The closedness means that, for 
any positive integer r and any r-index k, the time derivative /i^ can be expressed in terms of the mixed 
moments of order r and lower. This property is a corollary of the following general result. 
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Lemma 3: For any positive integer r and any r-index := (fci , . . . , fc^) G { 1 , . . . , ny, the mixed moment 
lAjc from (|38T ) for the plant observables, governed by the QSDE (fT9l) . satisfies 

= £ ^{^k,...kj_iFkj^kj+,...k,) 
7=1 

+ 52 H ^i'^ki...kj_igkjs'^kj+i...ke_igkfu'^ke+i...kr)- (39) 

Here, Fp is the />th entry of the drift vector F, and gps denotes the {p,s)th entry of the dispersion matrix 
G. 

Proof: By applying a multivariate version of the quantum Ito formula to and using the QSDE 
(fT9l) together with the quantum Ito product rules [13] (dt)^ = 0, {dt)dW = and ©, it follows that 

dZk = Y,^k,...kj_,dXkj 

7=1 

X (Sa,.^i...A:,+ SA,.^j...fe^_j(dX^^,)S^/'+i...-tJ 

= £ ^ki...kj_^ ( (Ffc.d? + ^fc..dW)E^.^j...fc^ 

m r , 

+ £ (Osugkjs £ '^kj+i...kf_igkfu'^kf+i...k,\ (40) 

where gp, denotes the pth row of the dispersion matrix G. The ODE (|39| ) can now be obtained by 
averaging both sides of (|40] ) and using the special structure ([5]) of the quantum state p(0). ■ 
It follows from Lemma |3] that if F is an affine function of X and the dispersion matrix G is constant, as 
in the case where the plant-field coupling operators are linear in X, satisfying the CCRs, then the right- 
hand side of the ODE (l39l) is a linear combination of the mixed moments jUy, where the multiindices 
V := (vi, . . . , v^) have dimensions d ^r. More precisely, 

r 

Afc=£ £ V^fc,vMv, ke{l,...,nY, (41) 

d=Ove{l,...,ny 

where i//yt,v complex numbers which are found from (l39T l. and the convention /i^ := 1 is used for the 
moment associated with the 0-index. Equivalently, an infinite dimensional vector /i := {lJ-k)ke{i ny.r^o^ 
formed by /i^ and the mixed moments /i^ for all possible n'' multiindices k E {I ...,nY of orders 
r — 1, 2, 3, . . ., satisfies a system of linear ODEs fl = ^'jl, where ^ := ( V^a:,v) is an infinite-dimensional 
block-lower triangular matrix. The diagonal block of ^ associated with the moments of order r is a 
matrix of order n''. Hence, the solution of the system of ODEs (|4T1) can be represented as /i(?) = 
e^'ll{0), provided all the moments of the initial plant state 07(0) are finite. The matrix exponential 
e*^' is practically computable due to the block-lower triangular structure of ^. Thus, the algebraic 
closedness (|4T1) allows the system of linear ODEs for the moments (l38l) to be integrated (numerically 
or analytically) recursively with respect to r, starting from the mean values of the plant observables for 
r = 1 . In particular, the mean vector and the quantum covariance matrix 

a:=EX, S:=co\{X)=E{XX^)-aa'^, (42) 

with the latter consisting of central moments of second order, satisfy the ODEs 

a=Aa + 07, S =AS + SA^ ^Bm'^ , (43) 
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which allow the steady-state values of the moments to be found from the appropriate algebraic equations 
if the matrix A, defined in (l37l) . is Hurwitz. Now, a similar reasoning shows that the moment dynamics 
(|39| ) retains the algebraic closedness (1411 in a more general case where both the drift vector F and 
the dispersion matrix G of the QSDE (fT9l) are qffine functions of X (in the linear case above, G was 
constant). This corresponds to a wider class of open quantum systems introduced in the next section. 

VI. QUASILINEAR OPEN QUANTUM SYSTEMS 

Retaining the assumption of Section |IV] that the plant observables satisfy the CCRs (l27l) . we will 
now consider a wider class of plant-field interactions in which the coupling operators are 
quadratic polynomials of the plant observables: 

hj = {Mj + yJ/2)X, Yj := RjX. (44) 

Here, Mj denotes the jth row of a matrix M E W^^'\ which describes the linear part of the coupling 
as in (l30l) . and G S„ are given matrices which specify the quadratic part. An equivalent 

vector-matrix form of (|44l i is 

h = (M + yT/2)X, Y:=[Yi ... y„] , (45) 

where Y is an {n x m) -matrix with columns , . . . , Y,n whose entries are linear combinations of the plant 
observables. In the case of quadratic plant-field coupling (l44l) . an argument, similar to the derivation of 
(|36l ) from (|35] ). allows the kth column (l22l) of the dispersion matrix G from (|20] ) to be computed as 

gk = i[{Mk + X^Rk/2)X,X] = @{Mj + Yk), (46) 

where MJ is the kth column of the matrix M^, which, in view of (|451) , implies that 

G = ®{M^ + Y)=B + ®[RiX ... R^X], (47) 

where the matrix B is defined by (l33l) . Therefore, the entries of G are affine functions of the plant 
observables. From (|44] ) and (|46] ). it follows that the contribution of the operators 

i[hj,g,] =i[iMj + Yj/2)X,mkX] 

=imk[{Mj+X^Rj/2)X, X] = @Rk@{Mj+Yj) (48) 

to the right-hand side of (1211 is linear with respect to the plant observables: 

m m 

i £ (Djk[hj,gk] = £ (Djk@Rk@{M] + RjX). (49) 

j,k=l j,k=l 

Thus, in the case of canonically commuting plant observables and quadratic plant-field coupling (|44] ). the 
dispersion matrix G is an affine function of X, while Sf{X), given by (|2TI) . is a cubic polynomial of X. 
The latter property suggests finding a Hamiltonian H in the form of a quartic (degree four) polynomial of 
the plant observables such that the corresponding cubic polynomial i[H,X] counterbalances the quadratic 
and cubic terms in =Sf(X), thus making the drift vector F in (|20l) an affine function of X: 

F=AX + I5, (50) 

where A G R"^" and /3 G M". Together with G depending affinely on X as described by (|47t . the resulting 
quantum plant, governed by the QSDE 

dX = {AX + I5)dt + @{m'^ + Y)dW 

m 

= {AX + p)dt + £ {Mj + RjX)dWj, (51) 

7=1 

which we will refer to as a quasilinear open quantum system, retains the algebraic closedness (|4ll in 
the moment dynamics (|39l ) as discussed in Section |Vl In the next section, we will show that the problem 
of finding such a Hamiltonian H is simplified significantly by taking physical realizability conditions 
into account. 
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VII. PRESERVATION OF CANONICAL COMMUTATION RELATIONS 

The commutator of two observables r/(?) and ^{t) on the product space M'®^ inherits the evolution 
^ with the unitary matrix U{t): 

\^{tu{t)\=v{t)\m.m\vit). (52) 

Hence, if r](0) and ^(0) satisfy a CCR, that is, if [t](0),^(0)] = (p is the identity operator /^®^ up 
to a complex multiplier 9, then the unitarity of U{t) and dSll) imply that \r\{t), C,{t)\ = (pU{t)^U{t) = 
[77(0), ^(0)] for all t^O. Therefore, any CCR between observables on the space J^®^ is preserved 
in time. In particular, if the plant observables Xi(0),. . .,X„(0) are in CCRs with each other, then, the 
preservation of these CCRs is a necessary condition for physical realizability (PR) of a QSDE of the 
form (fT9l ). Here, in accordance with [5], [12] in the linear case and [9] for nonlinear systems, PR is 
understood as existence of a plant- field energetics model, specified by the pair {H,h), which generates 
the particular drift vector F and dispersion matrix G as described by (|20l ). We will now obtain CCR 
preservation conditions for the quasilinear quantum plant governed by the QSDE (ISTI) which corresponds 
to (fT9l ) with the drift vector F and dispersion matrix G given by (|50l) and (l47l) . respectively. To this 
end, let ^ and M denote linear operators which map an n-dimensional vector u and a matrix T of 
order n to two matrices of order n: 

m 

J^{u):=@ £ JjkiMju^Rk + RjuMk)Q, (53) 

m 

M{T):=® £ JjkRjTRk®. (54) 

The significance of these operators is clarified by the following lemma which is instrumental to the 
proof of Theorem [T] 

Lemma 4: The matrix GJG^, associated with the dispersion matrix G in (|47l) . is a constant complex 
matrix (independent of X and the initial quantum state of the plant) if and only if the operators ^ and 
defined by (l53l) and (l54l) . both vanish on M" and S„, respectively. In this case, 

GJG'^ = BJB^ - /^(0)/2, (55) 

where the matrix B is defined by (l33l) . 

Proof: By combining ([33]), (gV]) with ([53]), dSH), it follows that 

G7G^ = - 0(M^ + Y)J{M + Y^)@ 

=BJB^ - 0(M^7y^ + YJM)@ - @YJY'^@ 
=BJB^ - J^iX) - M{XX"^) 

=BJB^ -^{X)-M{U)- jl. (56) 

where the operator matrix XX^ is split into a symmetric part 15 and the antisymmetric part [X,X^]/2 ~ 
i®/! as 

XX^ = O + /0/2, O := {XX^ + (XX'r)^) /2. (57) 

Hence, if ^ = on R" and =^ = on §„, then the terms ^{X) and M{iy) vanish in (|56l) . thus proving 
the sufficiency of these conditions, and the relation (|55T ) follows. In order to prove the necessity, suppose 
GJG^ is a constant complex matrix. Then it coincides with its expectation over any quantum state where 
X has finite second moments, that is, 

GJG^ = BJB^ -iM{®)/2-J^{a) -^(aa'^ + E), (58) 
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where a is the mean vector and 

E:=Re5 (59) 



is the real part of the quantum covariance matrix of X from (|42]) . and use is made of the relation 
= E + aoc^ which follows from (|57T ). By using a Gaussian initial quantum state [14] for X with 
an arbitrary mean vector a and the quantum covariance matrix 5 = E + /0/2, where E G §„ is varied 
independently of a subject to 5 )^ 0, it follows from the constancy of the matrix GJG^ in (l58l) that the 
linear operators ^ and vanish on the spaces R" and §„, respectively. This establishes the necessity 
and completes the proof. ■ 
Since is antisymmetric, the matrix ^(0) in (l56l) does not have to vanish under the assumption 
of Lemma |4] that =^ = on S„. A sufficient condition for this assumption to hold can be obtained by 
using the vectorization of matrices [10]: 



£ jjk{mj)®{mk)=Q. (60) 



More precisely, the condition (1601 ) is necessary and sufficient for the operator ^ in (|54l) to vanish on 
the space M"^" which contains S„. For what follows, we define, in a similar fashion to (|53] ). (|54|) . linear 
operators (S and Y which map an n-dimensional vector u and a matrix T of order n to two matrices of 
order n: 

m 

<§{u):=@ £ Vjk{M]u^Rk + RjuMk)@, (61) 

m 

r{T):=& £ vjkRjTRk®, (62) 

j,k=l 



where vjk are the entries of the real part of the quantum Ito matrix from 

Theorem 1: The quasilinear QSDE (ISTI) preserves the CCR matrix of the plant observables from 
(|27]) if and only if the linear operators ^ and M, defined by (l53l) . (|54|) . vanish on R" and §„, respectively, 
and 

A0 + + BJB^ = r (0) . (63) 
Proof: By combining the quantum Ito formula with the bilinearity of the commutator as d[X,X^] = 
[dZ,X^] + [X,dX^] + [dX,dX^], and using ([T9l ) together with the quantum Ito product rules, it follows 
that 



d[x,xT] = [Fdt + Gm,x^] + [X, {FAt + cm)^] 
+ [cm, [cmf] 

f m . 

= ([F,xT] + [X,fT]+/G7GT+ £ 

m 

+ £(b,^^] + [^,^fc])dWfc 

k=l 

m 

=i(A@ + @A^ + GJG^-@ £ w;^i?^0i?^0jdr. (64) 

j,k=l 



10 



Here, use has also been made of the commutativity between dW and the adapted processes F, G, X 
and the relationship 

m 

[GdW,{GdW)^]= £ [gjdWj,gJdWk] 

m 

= E [[Sj^sl]0ijk+i8j8pjk - i[gj:gJ]Jjk)dt 

j,k=l 

112 

=i(GJG'^~® £ coJj^Rj®Rk®]dt, (65) 

where Jji^ = llmcOjk is the ( j, k)th entry of the matrix J from (|3]), so that cOjk — iJjk = 0}jk- In turn, (|65] ) 
follows from the commutator identity [(pi//, ar] = [(p, a] i/ZT + ^[1//, a]T + a[(p, t] 1//+ C7(p[v/, t] which 
reduces to 

[(PV/, ctt] = [(p, a] v/T + C7(p[v/, t] 

= [(p,a]V/T+(pC7[V/,T]-[(p,(7][VA,T], (66) 

provided [i//, a] = and [9, t] = 0. More precisely, (|66] ) is applied to the case where ^, o are entries 
of the dispersion matrix G, while i//, T are those of dW. Note that the following term in (|64l) 



vanishes regardless of the particular form (I46l l of the columns of the dispersion matrix G in the case 
of quadratic plant-field coupling (|44|) . In fact, (l67l) follows from the general definition (|22] ) and the 
Jacobi identity [11] combined with the CCRs ^ whereby [X, [hk,X^]] + [[hk,X],X^] = [hk, [X,X^]] = 
i[hk,&] = 0. Also, in dMl), we have used the relation [gj.gj] = -&Rj[X,X'^]Rk& = -iQRj&Rk& which 
follows from (|46] ). Now, the CCRs (|T7l) are preserved in time, that is, the left-hand side of the QSDE 
(|64|) vanishes identically, if and only if so does the right-hand side: 



Ae + ®A^ + GJG'^ -® £ WJj^Rjmk@ = 0. (68) 

j,k=l 

The fulfillment of (|68l) is only possible if GJG^ is a constant complex matrix. By Lemma IH this 
property is equivalent to that ^ = on M" and ^ = on §„, in which case GJG^ is given by (|55] ). 
Now, since ®LJ,k^i WJj^RjQRk® = ^(0) - /2 in view of W~j^ = vjk - iJjk/2 and dH), then 

substitution of (ESI) into dMl) yields A0 + 0A'^ + 575'^ - ^(0) = 0, which is equivalent to (l63]), and the 
proof is complete. ■ 

Therefore, Theorem [T] imposes constraints (in terms of the operators ^ and M) which the quadratic 
plant-field coupling operators (|44]) have to satisfy in order to make an affine drift term of the QSDE 
(ISTI) achievable through an appropriate choice of the plant Hamiltonian H. 

Lemma 5: Suppose the CCR matrix in (ITT]) is nonsingular, and the operators ^ and ^ in (1531) . 
(|54l) . associated with the quadratic plant-field coupling model (|44l) . satisfy the conditions of Theorem [T] 
Then the GKSL vector ^{X) in (|2TI) is a quadratic polynomial of the plant observables with the leading 
term @YJMX/4, that is, 

^{X)= eYJMX I A + (affine function of X) . (69) 
Proof: Since det0 ^ 0, the relation h = (A/0 — G^)0^^X/2 between the vector h of quadratic 
plant-field coupling operators (|45T ) and the corresponding dispersion matrix G in (l47l) implies that 

GJh = GJ{M& - G^)&^^X/2 

= {@M^JM - GJG^&^ + &YJM)X/2, (70) 
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where use is also made of GJM = &{M +Y)JM. Now, if the operators ^ and ^ vanish on M" and 
S„, then, by Lemma HI the matrix GJG^ is constant. Note that for an arbitrary quadratic plant-field 
coupling model without PR constraints, GJG^ would be quadratic and GJh in (iTOl) would be a cubic 
polynomial. Therefore, the constancy of GJG^ reduces GJh to a quadratic polynomial of X, with its 
leading (quadratic) term being &YJMX/2 in view of the linear dependence of 7 on X. It now remains 
to substitute (TTOI) and (|49l) into (|2T1) in order to verify that ^{X) is a quadratic polynomial of X with 
the leading term QYJMX/4 as in (|69| ). where the calculation of the affine part is omitted for the sake 
of brevity. ■ 
Lemma [5] suggests that the class of candidate plant polynomials H for counterbalancing the nonlinear 
terms of the GKSL operator ^{X) by i[H,X] (to achieve an affine drift vector in the governing QSDE) 
can be reduced to cubic polynomials. One of such Hamiltonians is found in the next section. 

VIIL CUBIC PLANT HAMILTONIAN 

The following theorem provides a characterization of the class of quasilinear quantum stochastic 
plants described by the QSDE (ISTI) associated with the quadratic plant-field coupling model (|44l) . For 
its formulation, we introduce a Hamiltonian 

H :=Y^X+X'^RoX/2-X^YJMX/n, (71) 

which is a cubic polynomial of the plant observables. Here, yG M" and Rq E S„ are arbitrary vector and 
matrix which specify the quadratic part of H, while 

A :=X"^YJMX = [X^RiX . . . X^R,nX] J MX (12) 

is a homogeneous cubic polynomial of X, specified by the parameters M E M™^" and . . .,Rm E E>n 
of (gH)-®. 

Theorem 2: Suppose the observables of the open quantum plant under consideration have a non- 
singular CCR matrix in (|27l) . and the plant-field coupling is described by the quadratic model (|44] ) 
whose parameters satisfy the conditions of Theorem [T] Then the cubic plant Hamiltonian H, described 
by (ItTI) . is self-adjoint and leads to a quasilinear QSDE (|5TI) . 

Proof: Since the quadratic part of H, which was discussed in Section |IVl is a self-adjoint operator 
which contributes an affine function of X to the drift vector F in (|20|) . we will consider the cubic part 
of (TTTI) . In order to verify that A in (|72|) is indeed a self-adjoint operator, note that the definition of Y 
in (|44)) - (|45l) implies that the entries of the matrix YJM are linear combinations of the plant observables 

YJM=j^^^Xk, (73) 

k=l 

whose coefficients comprise matrices 4>i, . . . ,4>„ G M"^" as 

m 

4>,= £ Jp,{Rp).kM,. (74) 

Here, {Rp),k denotes the ^th column of Rp, and Ms is the 5th row of M as before. From the definition 
(|53l) of the operator it follows that, if det0 ^ 0, then the condition of Theorem [T] that ^ vanishes 
on M" is equivalent to the symmetry of all the matrices 4>yt in (1741) . Hence, substitution of (1731) into 
m yields At = {Zlk,i=ii^k)jiXjXkXey = L';,,,^=i(4>^),^XA^, = A, where (4>,),^ = (4>,)^, denotes 
the (j,^)th entry of <I>yt G S„. We will now compute the contribution of the cubic term A of H from (fTll 
to i[H,X]. To this end, by representing the operator A in (1721 ) as 

A = Z^JMX, Z:=Y^X = (X^RkX) (75) 
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where the m-dimensional vector Z consists of self-adjoint operators, it follows that 

i[A,X] = i{Z^JM[X,X"^])^ - i[X,Z^]JMX 
= -(Z^JM®)^ + 2@YJMX 

= -OM^yy^X + 2®YJMX = 3®YJMX. (76) 

Here, use is made of the CCRs (ITTI) , the relation — z[X,Z^] = 207 is obtained from (|75] ) by regarding 
the entries of Z as quadratic Hamiltonians, and the symmetry of the matrix YJM = —M^JY^ follows 
from (TTTt . Now, (TtTT ) and (1761 ) imply that the contribution of A to i[H,X] is described by — /[A,Z]/12 = 
—&YJMX/4, which is the negative of the leading quadratic term of the GKSL operator =Sf (X) computed 
in (|69l) of LemmalU Therefore, —i[A,X]/12 + ^{X) is an affine function of X and so is the drift operator 
F in (|20l) which corresponds to the cubic Hamiltonian (TtTT ). ■ 
Thus, Theorem |2] constructs a physically realizable quasilinear quantum stochastic plant from an 
appropriately constrained quadratic plant-field coupling model and the corresponding cubic plant 
Hamiltonian with an arbitrary quadratic part. The matrix A in (|5TI) . whose calculation is omitted for the 
sake of brevity, can be recovered from the proofs of Lemma |5] and Theorem |2l 

IX. QUADRATIC STOCHASTIC STABILITY 

Due to the algebraic closedness in the moment dynamics (|4TI) . which extends from the linear case of 
Section |W] to the quasilinear quantum systems (ISTI) . the stochastic stability of such systems is amenable 
to analysis at the level of moments of arbitrarily high order. We will discuss the quadratic stability 
which is concemed with the first two moments. 

Theorem 3: For the quasilinear quantum stochastic plant governed by the QSDE (ISTI) and satisfying 
the CCR preservation conditions of Theorem \T\ the mean vector a from (|42] ). and the real part E of 
the quantum covariance matrix in (|59] ) satisfy the ODEs 

a=Aa + /3, (77) 
E = AE + EA^ - r (E) + BVB'^ + /4 

-<^{a)-y{aa^), (78) 

where the linear operators £' and Y are defined by (|6T| ). (|62l ). 

Proof: Both (1771 ) and (1781 ) can be obtained by specializing the general moment dynamics (|39l ) 
to the quasilinear case of affine F and G. Alternatively, (1771 ) is established by averaging both sides of 
(ISTI) and using the special structure dS) of the quantum state p(0). In a similar vein, by averaging the 
quantum Ito differential d(XX^), it follows that the matrix 

n{t):^E{X{t)Xit)'') (79) 

of second moments of the plant observables satisfies the ODE 

n=AU + UA^ + l5a'^ + a[5'^ + E{GaG^). (80) 

The representation Q. = V + iJ/2 of the quantum Ito matrix Q. from © implies that 

Oao'^ =GVG^ + iGJG^/2 = BVB^ - c?(X) - 'tiXX'^) 

+ i{BJB'^ -iM{@)/2)/2, (81) 

where we have also used (IMT ). (|62l) . Lemma|4]and Theorem [T] The averaging of (ISTT ) yields E(GflG^) = 
BVB'^ + M{®)/A-S'{a)--r{n)+iBJB'^/2 whose substitution into ^ leads to 

fl =An + IM^ + [5a'^ + a[5'^ + BVB^ + ^(0) /4 

- ^(a) - r (n) + iBJB^/2. (82) 
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Since the matrix 11 from (|79l ) is representable as n = E + oca^ + z0/2, then ([78] ) is obtained by taking 
the real parts on both sides of (|82] ) and combining the result with (1771) . ■ 
Theorem |3] shows that, unlike the mean-covariance dynamics in the linear case (|43T ). the Hurwitz 
property of the matrix A is sufficient only for the stability of the quasilinear quantum plant (|5TI) at 
the level of the first order moments. In view of (TTSl ). such a plant is quadratically stable if the real 
parts of the eigenvalues of the linear operator E i— )■ AE + EA^ — "^(E), acting on the space §„, are all 
negative. In this case, the steady-state values limy^+oo and limf^+ooE(?) are unique solutions of the 
corresponding algebraic equations obtained by equating the right-hand sides of (TTTI ). (l78l) to zero. Also 
note that, for the quasilinear quantum plant, the mean vector a influences the evolution of E by entering 
the right-hand side of (|78T ) in a quadratic fashion, whereas the dynamics of the mean and covariances 
in the linear case (|43l) are completely decoupled. 
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